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We report extended simulation results and their new analysis in two important gauge theories with 
twelve fermion flavors in the fundamental SU(3) color representation and two fermions in the 
sextet representation. We probe the Nf = 12 model with respect to the conformal window using 
mass deformed finite size scaling (FSS) theory driven by the fermion mass anomalous dimension. 
Our results at fixed gauge coupling show problems with the conformal scenario of the Nf = 12 
model. In the sextet model with two flavors, under the conformal hypothesis, we determine large 
values for the anomalous fermion mass dimension with y > 1. Since our sextet analysis favors 
the chiral symmetry breaking hypothesis without conformality, the large exponent y could play 
an important role in understanding the composite Higgs mechanism. The new results discussed 
here include our extended data sets and exceed what was presented at the conference. 
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1. Introduction 

New physics at the Large Hadron Collider could be discovered in the form of some new 
strongly-interacting gauge theory with a composite Higgs mechanism, an idea which was outside 
experimental reach when it was first introduced as an attractive scenario beyond the Standard Model 
[1-9]. The original framework has been expanded by new explorations of the multi-dimensional 
theory space of nearly conformal gauge theories [10-18] where systematic and non-perturbative 
lattice studies play a very important role. Interesting models require the theory to be very close 
to, but below, the conformal window, with the gauge coupling slowly evolving over a large en- 
ergy range. The non-perturbative knowledge of the critical number of flavors NJ lt , separating the 
conformal phase from the phase of the composite Higgs mechanism with chiral symmetry break- 
ing (#SB), is essential and this has generated much interest with many old and new lattice stud- 
ies [19-60]. The position of the conformal window with respect to the much discussed model of 
twelve fermions in the fundamental representation remains controversial with recent efforts from 
several lattice groups [19-34]. The position of the Nf = 2 sextet model with respect to the confor- 
mal window also remains unsettled [20,44—48]. 

Probing the conformal and #SB hypotheses we use two different strategies to deal with finite 
volume dependence. The first strategy extrapolates the spectrum to infinite volume at fixed fermion 
mass m where the leading finite size correction is exponentially small and determined by the lowest 
mass which has pion quantum numbers. From the mass spectrum of the infinite volume extrapola- 
tion we can probe the mass deformed conformal scaling behavior and compare with %SB behavior 
when the fermion mass is varied in the infinite volume limit. The second strategy takes full advan- 
tage of the conformal FSS behavior without intrinsic scale when pressing against the m = critical 
surface at fixed finite size L. Different from the first strategy, the finite volume corrections are not 
exponentially small and a much larger data set is analyzed closer to the critical surface. This is used 
in the Nf = 12 model, significantly extending our previously reported results [21]. We will also 
briefly summarize our main results in the sextet model where the first strategy is sufficient, since 
only runs at the lowest fermion mass show consistent and detectable finite volume dependence. 

We have used the tree-level Symanzik-improved gauge action for all simulations in this paper. 
The conventional /3 =6/g 2 lattice gauge coupling is defined as the overall factor in front of the 
well-known terms of the Symanzik lattice action. Its value is /3 = 2.2 for all simulations reported 
here for the Nf = 12 model. In the sextet model results are reported at /3 = 3.2. The link variables 
in the staggered fermion matrix were exponentially smeared with two stout steps [61]; the precise 
definition of the staggered stout action was given in [62]. The RHMC and HMC algorithms were 
deployed in all runs. For molecular dynamics time evolution we applied multiple time scales [63] 
and the Omelyan integrator [64]. Our error analysis of hadron masses used correlated fitting with 
double jackknife procedure on the covariance matrices. The time histories of the fermion conden- 
sate, the plaquette, and correlators are used to monitor autocorrelation times in the simulations. 

2. Twelve fermion flavors in fundamental SU(3) color representation 

Extending our earlier work [21], we have new simulation results at j3 = 2.2 in the fermion mass 
range m = 0.002 - 0.025 at lattice volumes 20 3 x 40, 24 3 x 48, 28 3 x 56, 40 3 x 80, and 48 3 x 96. 
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The extended data base now spans the m = 0.002 — 0.035 range. The new lowest fermion mass 
runs at m = 0.002,0.004,0.006,0.008 can be used in the conformal FSS analysis which over the 
full set would correspond to a variation of the pion correlation length in the 2.5 to 20 range in the 
infinite volume limit. Results from the two lowest masses at m = 0.002,0.004 are not included in 
the current analysis and will be reported later. For further control on finite volume dependence, 
large 48 3 x 96 runs were continued to two thousand trajectories at m = 0.01 and m = 0.015. Four 
runs were further added at 40 3 x 80 with m = 0.01,0.15,0.02,0.025. The new and refreshed data 
set was subjected to conformal FSS analysis and ^SB tests of the (yy) chiral condensate. 

2.1 The phase diagram in the j8 — m plane 

The phase structure of the model remains controversial, particularly the critically important 
weak coupling phase. In addition to our spectroscopy and conformal FSS runs, we ran extensive 
scans at various fixed volumes and fixed fermion masses to explore the bulk phase structure. The 
bare coupling /3 was varied over a large range starting from very small j8 values deep in the strong 
coupling regime to the weak coupling phase at /3 = 2.2 where the conformal and #SB analyses 
were done. Fermion masses m = 0.007,0.01,0.02 were used in the scans with spatial lattice sizes 
L = 8, 12, 16,20,24,32 running a large densely spaced set in the important and much discussed 
intermediate region in transit from strong coupling to weak coupling. These scans were also ex- 
tended to Nf = 2,4, 6, 8, 10, 12, 14, 16 flavors. We will briefly summarize next what is known about 
the bulk lattice phase structure. 




P P 

Figure 1: On the left, scans of the phase diagram by monitoring the chiral condensate are plotted as a function of j8 at 
two different fermion masses. The schematic bulk phase diagram is sketched on the right. 

Two representative scans of the bulk behavior of the chiral condensate (yy) are shown in Fig- 
ure 1 as we vary j8 from strong to weak coupling. Three distinct regions emerge at fixed volume and 
fixed fermion mass showing strong coupling behavior for /3 < 1.4 with a large chiral condensate, 
an intermediate phase for 1.4 < j8 < 1.8 with sudden drop in (yip), and a weak coupling phase 
for /3 > 1.8 with further drop in (yy). A similar structure of three regimes was also seen in scans 
at Nf = 8. Our physics simulations were done well inside the weak coupling phase at j8 = 2.2 as 
indicated in Figure 1. A similar structure has been observed independently by Deuzemen et al. [28] 
and Cheng et al. [31]. The newfound order parameter of broken shift symmetry in the intermediate 
phase is the most interesting development in the study of the esoteric intermediate phase [31]. It 
only exists in a finite interval of the lattice gauge coupling for small enough fermion masses, as 
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schematically sketched in Figure 1. The real interest is, of course, in the nature of the weak cou- 
pling phase. Based on axial U(l) symmetry considerations, arguments were presented in [28] in 
favor of conformal symmetry in the weak coupling phase. This argument was criticized and refuted 
in [31] based on new details of the broken shift symmetry with chiral symmetry restoration they 
discovered at zero temperature in the bulk intermediate phase. 

Cheng et al. also presented their first weak coupling results on the Polyakov loop, the chiral 
condensate, and spectroscopy as indications of conformal symmetry in the weak coupling phase. 
The blocked Polyakov was reported to jump from zero to a large finite value in crossing to the 
weak coupling phase [31]. A confining potential was found in the intermediate phase with broken 
shift symmetry which turned into Coulomb potential without string tension in the weak coupling 
phase [31]. It was also noted that the observed chiral condensate and the related Dirac spectrum 
show the recovery of exact chiral symmetry in the massless fermion limit of the weak coupling 
phase consistent with observed degeneracies of parity partners even at finite fermion masses. 

The results in [31] suggesting a chirally symmetric deconfined conformal phase are in con- 
tradiction with what we reported earlier [21] and further confirmed in the extended new analysis. 
Using lattice volumes several times larger we find a vanishing Polyakov loop at zero temperature 
in the weak coupling phase and a confining potential at a pion mass which is lower than in [31]. 
We also find the parity partners split at finite fermion mass. Our findings in large volumes are 
consistent with a chirally broken weak coupling phase. As a first step to resolve the contradictions, 
large finite volume effects acknowledged in [31] will have to be brought under better control. 

In the next sub-sections we will briefly summarize our results on the chiral condensate, the 
finite temperature phase transition and tests of the conformal hypothesis in the weak coupling 
phase. 

2.2 Chiral condensate and #SB test 

In the extended new analysis the chiral condensate (Fig. 2) remains consistent with #SB in 
the massless fermion limit. Small changes in the fits are mostly driven by the two lowest fermion 
masses at m = 0.01 and m = 0.015 where runs on the largest 48 3 x 96 lattices were extended and 
new runs at higher masses on 40 3 x 80 lattices were added. A slight drift at the lowest m = 0.01 
fermion mass was detected in the connected part of the condensate even after 1,400 trajectories 
(i.e., MD time units) and the run was continued to 2,000 trajectories. Since the finite volume 
analysis is incomplete, the largest volumes are used in the fitting range of the fermion masses. 
Finite volume extrapolations have to be completed before definitive conclusions can be drawn to 
establish a non-vanishing condensate (yy) in the m — > limit. 

Details of the fitting procedure and the notation in Figure 2 were explained earlier [21]. The 
chiral condensate has a spectral representation [65] where the UV-divergent integral is written in 
a twice- subtracted form [66]. The UV contribution, which is divergent when the cutoff /x a 1 is 
removed, has a linear term « aT 2 ■ m and there is a third-order term «m 3 which is small and hard 
to detect for small m. The IR finite contributions to the chiral Lagrangian have a constant term 
mBF 2 , a linear term mB • m, a quadratic term ps B 3 F 2 ■ m 2 , and higher order terms, in addition to 
logarithmic corrections generated from chiral loops [67]. We kept a constant IR term and the linear 
and second order terms with UV and IR contributions. The second order fit in Figure 2 gives a non- 
vanishing condensate in the chiral limit which is roughly consistent with the GMOR [68] relation 
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Figure 2: The second order polynomial fit to the chiral condensate is shown on the left plot in subtracted form as 
explained in [21]. The middle plot is the quadratic fit to {yfyf) — m ■ % con directly measured from zero momentum sum 
rules and independently from functions of the inverse staggered fermion matrix. The right side plot shows the thermal 
history of the subtracted form of the condensate at the lowest fermion mass on the largest lattice. 

{xjfXjf) = 2F 2 B with the measured low value of F and the value of B from logarithmic fit to the 
Goldstone pion. The deficit between the two sides of the GMOR relation is sensitive to the fitting 
procedure and uncertainties in the determination of B. Trying to identify chiral logs is beyond the 
scope of our simulation range. For independent determination, we studied the subtracted chiral 
condensate operator defined with the help of the connected part % con of the chiral susceptibility % 
as defined in [21]. The removal of the derivative term significantly reduces the dominant linear 
part of the condensate. Although the two independent determinations give consistent non- 

vanishing results in the chiral limit, we cannot consider the results definitive. The drop of the chiral 
limit intercepts after extended runs is noted in comparison with earlier results [21]. 

2.3 Finite temperature transition 

We present some preliminary results from our extended studies of the finite temperature tran- 
sition. If the ground state of the model has #SB, a phase transition is expected at some finite 
temperature to restore the chiral symmetry in the limit of massless fermions. Based on universality 
arguments [69] the transition would be expected to be of first order. This is not entirely clear and 
warrants further investigations. On our largest 48 3 x Nt lattices, at fixed m = 0.01 and j3 = 2.2, 
the temperature was varied through an Nt sequence while the scatter plot of the Polyakov loop 
was monitored along the euclidean time (inverse temperature) direction in each run. The chiral 
condensate was a l so monitored in the runs. As the temperature is increased a clear sudden 

transition is observed in the Nt = 6 — 10 region where the Polyakov loop distribution jumps from 
the origin to a scatter plot with non- vanishing real part. It would be difficult to reconcile this jump, 
as shown in Figure 3, with conformal behavior in the zero temperature bulk phase. 

Although we have results for temperature scans at multiple gauge couplings, fermion masses, 
and spatial volumes, all consistent with a finite temperature transition, caution is necessary before 
firm conclusions can be reached. Confirming the existence of the #SB phase transition will require 
the m — > limits of (yy) an d the Polyakov loop distribution. The chiral condensate is a good 
order parameter for the transition. The Polyakov loop, like in QCD, could detect deconfmement 
in the transition with unsettled interpretation as order parameter. The behavior of the renormalized 
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Polyakov loop scatter plot at finite temperatures 
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Figure 3: The scatter plot on the left side shows the Polyakov loop in the time direction as the temperature is varied 
using a sequence of lattice sizes 48 3 x Nj at fixed m = 0.01 and j3 = 2.2 with Nj varied in the Nj = 6, 8, 10, 12, 16 range. 
On the right side plot the subtracted chiral condensate is plotted at j3 = 2.2 for a sequence of lattice sizes 48 3 x N T 

as 

the fermion mass is varied. The subtracted linear term 12.5 • m, chosen for the presentation of the data, is not determined 
from fits. 



Polyakov loop is consistent with the scatter plot of Figure 3. 

2.4 Conformal finite size scaling analysis 

The expected leading FSS form for any mass M, or for F K , scaled with the linear size L of the 
spatial volume, is given by a scaling function L-M = f(x) where x = L ■ m 1 ' l+ ? is the conformal 
scaling variable. The scaling form sets in close to the critical surface for small m values. The scaling 
functions fix) can depend on the quantum numbers of the states but the scaling variable is expected 
to have the same form with identical y exponent in each quantum number channel [39-42] . In sub- 
leading order there are conformal FSS scaling violation effects which are exhibited as a combined 
cutoff and L-dependent leading correction with the modified form L-M = f(x) +L~ w g(x) where 
the scaling correction exponent ft) is determined at the infrared fixed point (IRFP) g* of the /3- 
function as ft) = j3'(g*). This assumes that the mass deformation away from the critical surface is 
the only relevant perturbation around the IRFP. The leading scaling correction term close enough to 
the critical surface dominates any other corrections which are supressed by further inverse powers 
of L. To detect the leading scaling violation effect requires high precision data with fits to scaling 
functions f(x) and g(x) and the critical exponent ft). 

We applied conformal FSS theory to our data sets in the fermion mass range m = 0.006 — 0.035 
with lattice sizes ranging in the fits from 20 3 x 40 to 48 3 x 96. Two different FSS fitting procedures 
were applied. In the first procedure, we defined a scaling function fix) for each mass M with five 
independent fitting parameters. The fitting function fix) was divided into two regions separated 
at the joint x = x cut . Different forms were chosen on the two sides of x cut from the expected 
conformal behavior. For large x > x cut , the function f(x) = c^x + c exp {c\x)~ l l 2 e,x<p{~ c i x ) with 
parameters c\ and c exp describes the L-independent limit M ~ cim 1 / 1+1 ' at fixed m and L — > °°. 
The amplitude sets the size of the leading small exponential correction from the finite volume 
effect of the lightest Goldstone pion state wrapping around the spatial volume. Since /(0) = co 
is expected from conformal FSS with some power corrections at small x, we applied the simple 
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conformal FSS (Goldstone pion - PCAC channel) 



conformal FSS (Fpi - PCAC channel) 
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Figure 4: Conformal FSS fits in four different quantum number channels. The fits are performed in each 
channel separately. Since the y values vary considerably from channel to channel, a simultaneous global fit 
to the combined channels with the same / exponent, as required by conformal FSS theory, is bound to fail. 



ansatz f(x) = cq + c a x a for x < x cut (a more general polynomial function in the small x region 
is not expected to change the conclusions from the fits). From the fit to the PCAC Goldstone 
pion channel the parameter c n = c\ was determined and used as input in the exponential terms of 
the other channels with exp(— c^L). The critical exponent y was included among the five fitting 
parameters, in addition to Co, c\ , c exp , and x cut . 

The composite particle masses in several quantum number channels can be reasonably fitted 
with conformal scaling functions f(x) as shown in Figure 4 but the values of the critical exponent 
y are incompatible across different channels. The required global conformal FSS fit will fail with 
a single exponent y across all quantum numbers. In the fits for F n in the PCAC pion channel 
we only kept four parameters because the asymptotic form with exponentially small correction 
was zero within error. Actually, the data of F K did not allow a successful conformal fit with any 
shape chosen for its scaling function f{x) which looks very different from the scaling functions of 
composite particle masses. The unexpectedly curious behavior of the F n data set against conformal 
FSS remains unresolved. 
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2.5 Generalized FSS fitting procedure with spline based general B-form 

Following a new fitting strategy, we investigated if the failed global conformal FSS analysis 
can be attributed to restrictions on the conformal scaling functions f(x). The restrictions were 
manifest in the physics-motivated fitting procedure we applied above. Our new general approach 
is different from [24,45] but addresses related issues. We developed a general least-squares fit- 
ting procedure to the scaling functions using the B-form of spline functions without any further 
restrictions. In this procedure, the function f(x) is described by piece- wise polynomial forms con- 
structed from spline base functions with general coefficients in overlapping intervals of the scaling 
variable x. The shape of the B-form can be changed without limitations by increasing the number 
of base functions and the number of scaling intervals in x bracketing the overlapping data range. 
The details of this new analysis will be reported elsewehere [70]. 



conformal FSS fit with spline based B-form (Goldstone pion) conformal FSS fit with spline based B-form (cRho2 channel) conformal FSS fit with spline based B-form (F in PCAC channel) 




x= Lm* m x= Lmv m x= Lm* m 



Figure 5: Conformal FSS fits using spline based B-forms in three different channels. The fits are preformed in each 
channel separately with the question mark on y indicating difficulties of error estimates in bad fits of F n . 

Our fitting procedure in its setup requires two steps. In the first step, for any fixed choice 
of the exponent y, the best fitted function f(x) is determined in spline function B-form from the 
minimization of the weighted x 2 expression. According to a general algorithm, the x-range of the 
data set is divided into intervals separated by internal knots and adding end point knots for B-form 
spline construction. The number of coefficients is determined by the number of knots and the order 
of the spline polynomials of the sub-intervals. The weighted % 2 sum is minimized with respect 
to the coefficients of the base functions in the B-form. This will produce the best fit for fixed y 
with a minimized x 2 sum which will depend on y. In the second step, we minimize the % 2 sum 
with respect to y to determine the best fit of the critical exponent. The one-a confidence interval is 
determined from the variation of the % 2 sum as a function of y. 

In Figure 5 we show three typical fits for illustration. The fit to the Goldstone pion in the 
PCAC channel improved as expected, with considerable increase in the error. The tension across 
channels decreased, as illustrated by comparison with the rho-channel fit, but the fit to F n remained 
unacceptable. With the extended data set we are unable to reproduce results in [24,45] which used 
tables from our earlier limited subset of data [21] in favor of consistency with the conformal phase. 
It is important to emphasize that we have not reached definitive conclusions about the failure of 
conformal tests. As we stated earlier [21], we have not analyzed yet the leading scaling violation 
effects and did not investigate if the good scaling form in separate quantum number channels can be 
explained in the chirally broken phase by strongly sqeezed wave function effects. In disagreement 
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with [45], conformal FSS based analysis of the spectrum and related sum rules on moments of the 
correlators we have been developing are deep renormalization group based probes of the conformal 
phase. As explained in our forthcoming publication [70], we remain skeptical about the fitting 
procedure followed in [24] with efforts to rescue the conformal interpretation. The issues are not 
settled and ultimately will be decided in more definitive analyses. 

3. Two fermion flavors in the sextet SU(3) color representation 

This model has been studied recently by three BSM groups [20,44,46,48]. Our findings 
are different from results based on the Schrodinger functional [44, 46] and compatible with the 
finite temperature phase transition in [48]. The disagreement with Schrodinger functional results is 
particularly significant based on the lower bound y > 1 we find adopting the conformal hypothesis. 
This can be important in BSM applications and remains in contrast with the small exponent y < 
0.45 published in [46]. 

We have new simulation results at jS = 3.2 in the fermion mass range m = 0.003 — 0.010 on 
24 3 x 48 and 32 3 x 64 lattices. Five fermion masses at m = 0.003,0.004,0.005,0.006,0.008 are 
used in most fits. For further checks on finite volume dependence, a very large and expensive 
48 3 x 96 run was added recently at m = 0.003 to follow the strategy of finite volume extrapolation 
at fixed fermion mass m. We also have new preliminary simulation results at j8 = 3.25 in the mass 
range m = 0.004 — 0.008 on 24 3 x 48 and 32 3 x 64 lattices. Based on the chiral and conformal 
analyses of the model, continued runs at existing run parameters and new runs are planned at both 
couplings to further probe the conformal FSS hypothesis in the sextet model following the strategy 
we presented for the Nf = 12 model. 

3.1 Finite volume extrapolation 




Figure 6: Finite volume dependence at the lowest fermion mass for j8 = 3.2. The form of g\(\,r\) is a complicated 
infinite sum which contains Bessel functions and requires numerical evaluation [71]. Since we are not in the chiral log 
regime, the prefactor of the g\ (A , tj ) function was replaced by a fitted coefficient. The leading term of the function 
g\ (A, r/) is a special exponential Bessel function K\ (A) which dominates in the simulation range. 

Based on the #SB hypothesis, infinite volume extrapolations of the Goldstone pion, F n , and 
(\j7y) are shown in Figure 6 where g\ (A , 7] ) describes the finite volume corrections with X = M n ■ L 
and aspect ratio r\ = T /L from the lightest pion [72]. The fitting procedure approximates the 
leading treatment of the pion which wraps around the finite volume, whether in chiral perturbation 
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theory, or in Liischer's non-perturbative finite volume analysis [73]. This does not require to reach 
the 1-loop chiral log limit as long as the pion is the lightest state dominating the finite volume 
corrections. The infinite volume limits of M % , F„, and (\j/y) for m = 0.003 at /3 = 3.2 were 
determined self-consistently from the fitting procedure. Similar fits were applied to other composite 
states. Based on the fits at m = 0.003, one percent accuracy of the infinite volume limit is reached 
at M n L = 5. In the fermion mass range m > 0.004 the condition M K L > 5 is reached at L = 32. 
Although it will require high precision runs to test, we do not expect more than one percent residual 
volume dependence in the 32 3 x 64 runs for m > 0.004. Based on these observations, we will 
present chiral and conformal analyses with extrapolated infinite volume scaling behavior from the 
32 3 x 64 runs for m > 0.004. 

3.2 The chiral condensate and xSB 

We follow the analysis of the chiral condensate as described for the Nf = 12 model. The ( yy) 
condensate data were fitted with a third order polynomial of the form Co + c\m + C2m i while the 
condensate with derivative subtraction was fitted without the linear term. Both independently mea- 
sured quantities have to converge to the same chiral limit. The chiral condensate and its subtracted 



Nf = 2 sextet /3 = 3.2 (**) subtracted N, = 2 sextet (5 = 3.2 




m q 



Figure 7: For any given m > 0.004 the largest volume condensate data is used since the finite volume analysis remains 
incomplete. The two plots are discussed in the text. 

derivative version are shown in Figure 7 with a consistent strong %SB signal in the chiral limit. 

3.3 Spectrum and the #SB hypothesis 

Based on the analytic fermion mass dependence of the chiral Lagrangian, and using the lowest 
fermion masses in the m = 0.003 — 0.008 range, good polynomial fits were obtained without loga- 
rithmic loop corrections as shown in Figure 8 for four select states. The plotted 24 3 x 48 data points 
for m > 0.004 agree with the fitted data from the 32 3 x 64 runs indicating the infinite volume limit 
within the accuracy of the data. Small corrections, if required, should not effect the conclusions. 
Although we could fit M K and F % with the continuum chiral logarithms included, the separate sets 
of F and B from the fits are not quite self-consistent. A combined staggered SU(2) chiral perturba- 
tion theory fit is required for simultaneous fits of M % and F K with a consistent pair of fundamental 
chiral parameters F and B. The explicit cutoff dependent corrections to the F and B parameters 
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Figure 8: Polynomial fits from the analytic mass dependence of the chiral Lagrangian without logarithmic loop correc- 
tions are shown for the Goldstone pion, F x , M p , and the ++ state with mass Mjjiggs- The dashed line in the Goldstone 
pion plot shows the leading linear contribution. Fq on the top right plot sets the eletroweak vev scale. The disconnected 
diagram, which can shift the final value of the Higgs mass from Mq = 6.06 • on the bottom right plot, is not included 
in the calculation. 



would require further testing at weaker gauge couplings and using partially quenched staggered 
chiral perturbation theory. Our runs at /3 = 3.25 should provide the data for this analysis. 

3.4 Conformal hypothesis and the critical exponent y 

It is important to compare the polynomial fits with conformal scaling behavior for small mass 
deformations m. In the infinite volume limit the masses of composite particles and F K are expected 
to scale as M ~ with the same exponent y in all channels. When the four lowest fermion 

mass values closest to the critical surface are fitted separately with the leading conformal form, 
we get good % 2 fits but very different y exponents, which is not consistent with mass deformed 
conformal behavior. The conflicting fits are illustrated side by side in Figure 9 for the Goldstone 
pion and the F K decay constant. Fitting to the pion mass requires y = 1.091(34) while the F % fit is 
forcing y = 2.13(18). In the combined fit they compromise with y = 1.47(26) and the unacceptable 
X 2 /dof = 31.1. It is very difficult to see how this conflict, also in disagreement with [46], could 
be resolved within the conformal hypothesis. From the tests we were able to perform, the sextet 
model is consistent with %SB and inconsistent with conformal symmetry. It will require further 
investigations to show that subleading effects cannot alter this conclusion. We will consider com- 
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Figure 9: The left side plot and the middle plot represent separate conformal fits. The right side plot display the M x 
residuals from the global fit. It is unacceptable for F n as well. The global fit is trying to choose a y value between y ~ 1 
in the Goldstone pion channel and y ~ 2 in the F K fit resulting in a very high x 2 value. All fits are at j8 = 3.2. 

prehensive conformal FSS tests which do not rely on infinite volume extrapolation in the scaling 
fits. This is at a preliminary stage requiring new runs and systematic analysis. 

If %SB of the sextet model is further confirmed in the massless fermion limit, its relevance for 
the realization of the composite Higgs mechanism is transparent. The large anomalous exponent 
7 of our conformal fits will be interpreted in this case as an important ingredient of the model in 
the %SB phase. Importantly, the model has the perfect match of three Goldstone pions to provide 
the longitudinal components of the and Z bosons. To understand the slowly changing gauge 
coupling close to the conformal window without infrared fixed point will require high precision 
methods to calculate the renormalized gauge coupling and its beta function. This will demand 
extended and more reliable Schrodinger functional analysis or alternate methods which are being 
developed. The difference between the large exponent y reported here and the low value of y 
published earlier [46] is significant and will require clarifications. Conformal FSS tests very close 
to the critical surface will provide further independent checks of our results. 
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